Reconciling the efficacy and effectiveness of masking on epidemic outcomes

During the COVID-19 pandemic, mask wearing in public settings has been a key control measure. However, the reported effectiveness of masking has been much lower than laboratory measures of efficacy, leading to doubts on the utility of masking. Here, we develop an agent-based model that comprehensively accounts for individual masking behaviours and infectious disease dynamics, and test the impact of masking on epidemic outcomes. Using realistic inputs of mask efficacy and contact data at the individual level, the model reproduces the lower effectiveness as reported in randomized controlled trials. Model results demonstrate that transmission within households, where masks are rarely used, can substantially lower effectiveness, and reveal the interaction of nonlinear epidemic dynamics, control measures and potential measurement biases. Overall, model results show that, at the individual level, consistent masking can reduce the risk of first infection and, over time, reduce the frequency of repeated infection. At the population level, masking can provide direct protection to mask wearers, as well as indirect protection to non-wearers, collectively reducing epidemic intensity. These findings suggest it is prudent for individuals to use masks during an epidemic, and for policymakers to recognize the less-than-ideal effectiveness of masking when devising public health interventions.


Introduction
During the coronavirus disease 2019 (COVID-19) pandemic, universal mask wearing in public settings has been implemented in many places throughout the world as a means of mitigating the pandemic.A wide variety of masks, ranging from cotton masks to N95 respirators, have been used by the general population throughout the pandemic.Intuitively, masks create a physical barrier between the wearer's airway and the environment, reducing the chance of both disseminating and acquiring infection.In addition, a number of laboratory experiments, often based on manikin simulators, have demonstrated high efficacies of different masks including cotton masks; measured reductions of inhaled particles in laboratory settings ranged from approximately 50% to greater than 90% [1][2][3].However, the reported effectiveness of masking, demonstrated from randomized controlled trials (RCTs), has been much lower than the laboratory experiments, ranging from 10 to 25% [4,5] to inconclusive in a 2023 Cochrane review [6].Further, many places have endured multiple COVID-19 pandemic waves since March 2020, despite mask mandates implemented in public settings.The large discrepancy between the efficacy (measured under ideal conditions of mask quality and correct use, e.g. in laboratories measured by the filtration rate) and effectiveness (measured under real-world, often imperfect conditions, e.g. in RCTs) of masking has prompted debates on the need for mask use in the general population [7,8].
To explain this large discrepancy, studies have cited a number of factors.These mainly include compliance (e.g. the percentage of people actually wearing masks could be lower than the expected 100% among the test arm of an RCT), consistency (e.g.participants may not wear masks all the time, including the moment of exposure), and timing of implementation (e.g.there may be little circulation/presence of pathogen during an RCT to support demonstration of an effect) [8][9][10].However, it is difficult to test these factors in real-world settings and, as such, opinions and practices are divided-some dismiss the low effectiveness estimates from RCTs and continue to mask; others dismiss the high laboratory efficacies and refuse to mask.In addition, it remains unclear when and to what extent masking can be effective at the population level.
Here, to inform more effective individual preventive care and public health measures, we develop an agent-based model that comprehensively accounts for individual masking behaviours and efficacies as well as the nonlinear transmission process, and test the combined impact of these processes at the population level.Using realistic inputs of mask efficacy (i.e.50-80%, mean = 66%; electronic supplementary material, table S1) and contact data at the individual level, the model is able to reproduce the lower effectiveness (e.g. less than 25%) at the population level, as observed in real-world settings.Model results highlight the main factors contributing to the lower effectiveness, illustrate the complexity of measuring the underlying impact of masking due to nonlinear epidemic dynamics, and provide a more complete picture of the protective impact of masking (e.g.reducing repeated infections).

Overview of the agent-based model and analyses
The model used here has several key features.First, it explicitly simulates the transmission process to capture the contagious nature of infectious diseases.That is, individual risk of infection is determined by behaviour (e.g.contact rate and mask preference) as well as the people with whom an individual interacts (e.g. the infection status of contacts).This structure depicts the connectivity and spread of infection across a population, which is probably a more realistic representation of the real world (versus e.g.regression models).Second, the model separates household and non-household contact and transmission.This formulation not only allows modelling of the relative importance of household and non-household transmission based on contact data, but also allows modelling of specific masking behaviour in each setting (e.g.masks are often not used within a household).Third, the model records the history of each infection and vaccination to account for immune protection and waning for each individual and, through the aforementioned inter-person interactions, the combined outcomes at the population level.Fourth, the model is age-specific, which allows age-specific contact rates (e.g. higher among school-age children), efficacy of masking (e.g. higher for older adults) and disease severity (electronic supplementary material, table S1).This differentiation allows more realistic representation of population heterogeneity.Fifth, for a given mask type with known filtration efficiency, wearers may have a higher chance of infection when encountering a pathogen with a lower median infectious dose (ID 50 ) than that with a higher ID 50 , i.e. akin to having a lower mask efficacy.While we do not explicitly model pathogen-specific characteristics, the wide range of masking efficacies explored in our model simulations (electronic supplementary material, table S1) allows examination of infection outcomes (i.e.mask effectiveness) when mask efficacy is varied due to e.g.different IDs.Lastly, the model is run continuously over multiple years, where an individual is subject to varying background (environmental) infection risks as the population-level infection prevalence varies over time.As such, it represents mask effectiveness and the cumulative impact for different background (environmental) infection risks.Taken together, unlike previous approaches, we more comprehensively capture individual-level effects, nonlinear epidemic dynamics, and overall impact at the population level.
Using the model, we simulated a COVID-19-like pandemic and subsequent epidemic waves in a population of 100 000 people, from 1 March 2020 (roughly the pandemic onset) to 31 December 2024 (roughly a 5-year period).The basic reproduction number (R 0 , i.e. the average number of secondary infections in a naive population) was set to 2.6-similar to estimates for the ancestral SARS-CoV-2 [11][12][13][14] that initiated the COVID-19 pandemic-in our main analysis and counterfactual analysis using contact data from the POLYMOD study [15] and assuming an infection probability of 50% (10%) for physical (nonphysical) contact.To simulate pandemic dynamics, we used mobility data to adjust contact rates and vaccination data to account for vaccine protection, for all simulations.To test the impact of masking, we varied the masking rate (i.e. the fraction of population using masks, 10-100%) and the timing of masking (always wear masks outside the home versus when prevalence is above 0.01% to 2%) and compared simulated epidemic outcomes with a baseline setting without masking.
We compute the masking effectiveness at the population level as the reduction in a given health outcome (e.g. the total number of infections in a population) under a given masking setting compared with the baseline setting without masking (see Methods).In doing so, the population-specific characteristics are largely cancelled out, such that the effectiveness estimates would be more generalizable to different populations.This measure is also akin to those reported in cluster RCTs, which typically compare the overall attack rates in communities assigned to wear masks (where there may be different levels of adherence, similar to the simulations here under different masking settings) with communities not assigned to wear masks (similar to our baseline setting without masking).This measure is also consistent with non-cluster RCTs reporting reductions in the attack rate of a test arm with masking compared with a control arm without masking.In addition, as noted in the Introduction, reported masking effectiveness measures from various RCTs fall within similar ranges [4][5][6].As such, here we compare our model results with reported estimates in RCTs overall (i.e.without separating individual types of RCTs).

The model recapitulates masking effectiveness reported in randomized controlled trials
Figure 1 shows example simulated daily infection rates during the approximately 5-year study period.Under the most vigilant scenario in which people choosing to wear masks (i.e.wearers) always do so outside the home, increasing the masking rate greatly reduces epidemic size; this reduction is consistent across all years (figure 1a).However, even with a near 100% masking rate (i.e.all royalsocietypublishing.org/journal/rsif J. R. Soc.Interface 21: 20230666 except for those aged less than 1 year use masks), epidemic waves, albeit of much smaller size, persist.In addition, the effectiveness decreases substantially, when masking is implemented only when population prevalence of infection is high (e.g.greater than 1% in figure 1b).It is further evident from figure 2a, across a wide range of scenarios and time frames, that there are many instances that estimated effectiveness (measured as the reduction in total number of infections compared with the baseline without masking) is less than 25%-the higher bound of estimates reported in RCTs-including instances with high masking rates (e.g.greater than 80%) if masking is enacted late (see cells within the black polygons in figure 2a).

Factors contributing to lower masking effectiveness in the real world
To understand the main factors behind low masking effectiveness, we hypothesize that transmission may persist in household settings where contact tends to be greater, yet masking is rare.As shown in figure 3a (contact rates formulated using the POLY-MOD survey data [15]), for most individuals, household contact accounts for the majority of total contact, and are hence a likely principal source of infection.Thus, to test the impact of household transmission, we simulated a counterfactual setting in which people also wear masks at home (figure 3b).
Three main findings emerge from the counterfactual analysis.First, compared with when masking is solely applied outside the home (main setting, figure 2a), the effectiveness would be much greater should masking be also universally applied at home (counterfactual, figure 3b).For instance, assuming a 50% masking rate and masking at all times, over the 5-year study period, the counterfactual effectiveness is 62%, more than twice of the estimate in the main setting (25%; figure 2a, last panel, the cell in the first column and middle row).Second, it is evident from figure 3b that the effectiveness is higher than the expected direct effect should masking only affect the wearers.For instance, with a 50% masking rate and masking at all times, the expected effectiveness is 33% (i.e.50% × 66% mean efficacy), much lower than the 62% estimate.This greater effectiveness highlights the indirect protection of masking to the entire population including non-wearers, an effect that is hidden in real-world settings.Third, there is often a delay in awareness of disease prevalence and, in turn, a delay in masking.Such delays further facilitate introduction of infection into households and onward transmission (both within and outside households), leading to decreases in the overall Example simulated epidemic infection rates.For these simulations, the model assumes an R 0 of 2.6 and no mask wearing at home.Coloured lines (median) and shaded areas (interquartile ranges) show daily infection rates for simulated populations with different masking rates (see legend) and timing of masking: (a) mask wearers do so outside the home at all times or (b) only when infection prevalence is greater than 1%.For context, (c) shows changes in daily population mobility (representing reductions in contact rate due to social distancing rules) and seasonal infection risk (representing effects due to environmental conditions), used in the simulations.Mobility data for New York City during March 2020-June 2021 (shown here) were used directly as model inputs.From July 2021 onward, we assumed there was no social distancing and increased the mobility linearly to the pre-COVID level over a three-month period and then remained at that level (i.e.set to 1).Trends in seasonal infection risk were represented using a sinusoidal function (see details in Methods).
effectiveness over time.For instance, assuming a 50% masking rate and masking at all times, the estimated effectiveness in the main setting decreased from 44% during March-June 2020 to 25% at the end of the 5-year study period (figure 2a, first versus last panel), whereas the counterfactual estimates remained above 60% over the entire 5-year study period (figure 3b, first versus last panel).
The model results also illustrate how system nonlinear dynamics can affect masking effectiveness.For example, during the first wave (March-June 2020), estimated effectiveness surpasses the expected, under most combinations of masking rate and timing (figure 2a, the number in the cell > masking rate × 66% mean efficacy).This higher-than-expected effectiveness comes from the transmission reduction due to masking, which, in combination with the reduction in contact rate and less conducive environmental conditions in late spring/summer (figure 1c), reduces the effective reproduction number R e to less than 1, preventing a large pandemic wave (note: R e measures the average number of secondary infections and unity is the epidemic threshold).Conversely, when R e is much higher than 1, the effectiveness tends to be lower than expected due to 'leakage' via household transmission, as noted above.For instance, a large wave occurred during the 2020-2021 season (see figure 1a, second wave) under all masking rates, probably due to high population susceptibility and relaxed implementation of other non-pharmaceutical interventions (NPIs; e.g.social distancing; figure 1c, blue line) coinciding with more conducive transmission conditions during wintertime (figure 1c, red line).As a result, the reduction in epidemic magnitude was much smaller than during the first wave.This is also evident in the cumulative effectiveness over time (figure 2a; compare the higher effectiveness estimates for March-June 2020 in the first panel with the lower effectiveness estimates for subsequent seasons in remaining panels).

Measurement issues can further affect reported masking effectiveness
In addition to the factors noted above, the method and timing of measurement could both affect the reported effectiveness of masking.For the former, reported effectiveness would be lower than actual, if repeated infections are not counted.For instance, the number of people ever infected will inevitably increase over time.Studies that only measure the ever-infected while omitting    royalsocietypublishing.org/journal/rsif J. R. Soc.Interface 21: 20230666 reinfections will thus underestimate the actual (see the last four panels in figure 2b that discount reinfections, versus those in figure 2a that include reinfections).For the latter, prior epidemic dynamics could affect outcomes measured during a given period and consequently obscure the true effectiveness for masking.As shown in figure 1, epidemic size could vary across seasons, due to e.g.accumulated population immunity.Studies that measure the effectiveness based on a specific time period without accounting for prior dynamics could thus obtain vastly different estimates (figure 2c).These model results demonstrate potential measurement errors and highlight the importance of adopting more comprehensive measurements (e.g.including reinfections) and accounting for prior epidemic dynamics (e.g.examining cumulative effects).

Additional factors (sensitivity analyses)
The main analysis above assumes that infection risk via non-physical contact is 20% of that via physical contact (i.e.relative risk = 0.2).Depending on the etiologic agent, this relative risk could be higher or lower.In particular, higher infection risk via non-physical contact may apply to pathogens that transmit more efficiently via the airborne route.Thus, we tested the model under different nonphysical/physical contact-mediated relative infection risks (0.1, 0.5 and 1, versus 0.2 in the main analysis).Per the POLYMOD contact data, the basic reproduction number R 0 is slightly higher for a relative risk of 1 (3.2 versus 2.5-2.6 for other settings).As expected, it is more challenging to use masking to reduce infection when non-physical contact is as likely as physical contact to cause an infection (electronic supplementary material, figure S1, first row same risk versus other rows with lower relative risks).Nonetheless, the general findings noted in the main analysis are consistent with all other relative risk settings tested (electronic supplementary material, figure S1).
In addition, we hypothesize that mask wearers may be more risk adverse and thus may adopt additional precautionary measures such as reducing contact with others.We tested this by reducing the contact rate of mask wearers by 25%.As expected, the model estimates a further reduction in infections when mask wearers simultaneously reduce their contact rates (electronic supplementary material, figure S2 versus figure 2).For instance, the model estimates a 29% (versus 23% in the main setting) reduction in total infections over 5 years, if half of the population use masks when prevalence is greater than 0.01.
Lastly, over the course of the COVID-19 pandemic, viral mutations have led to more infectious SARS-CoV-2 variants.For instance, the Omicron BA.1 variant, which emerged in the end of 2021, had an R 0 around twice that of the ancestral SARS-CoV-2 variant [16].To examine the effectiveness of masking against more infectious pathogens, we tested the model under an R 0 of 5.2 (i.e. 2 × the value used in the main analysis, and similar to estimates for the Omicron BA.1 variant [16]).As shown in electronic supplementary material, figures S3 and S4, estimated effectiveness of masking decreases substantially under the much higher R 0 (e.g.estimated reductions of cumulative number of infections over the 5-year study period were half that estimated in the main analysis; see electronic supplementary material, figure S3a versus figure 2a).Similar to the first sensitivity analysis, these results are expected given the increased challenges to curb large epidemic waves caused by more infectious pathogens.Nonetheless, we note that reductions in infection for simulated populations wearing masks are still evident (electronic supplementary material, figures S3 and S4).     ) for the corresponding masking rate ( y-axis) and timing of masking (x-axis, percentages show population infection prevalence when masking is enacted).Bold fonts indicate the 2.5th% of the estimate is above 0 (i.e.significant at a = 0.05 level); note the 95% intervals here are probably narrower than measures in RCTs as the model assumes all infections are identified for the entire population.Negative values (i.e.increases; see grey cells) are possible due to model stochasticity (N.B., none are significant at a = 0.05 level).The black lines delineate cells with a median effectiveness estimate less than 25%.

Estimated effectiveness of masking at population and individual levels
To examine the effectiveness of masking at the population level, we compare simulated total numbers of infections in a population with different masking settings, with one with no masking.As shown in figure 2a and electronic supplementary material, table S2, at the population level, the effectiveness of masking depends on both the masking rate and timeliness of implementation.When implemented early during the epidemic, masking can substantially reduce infections (e.g.approximately 30% estimated reduction over the 5-year study period for a 70% masking rate in the main analysis assuming an R 0 of 2.6), despite the 'leakages' noted above.Model results using parameters based on COVID-19 also show similar reductions in hospitalizations (electronic supplementary material, figure S5 and table S2), suggesting similar downstream impacts on severe disease outcomes.However, timely implementation is crucial and there appears to be a threshold, after which estimated effectiveness declines rapidly-e.g. when community prevalence is greater than 0.2% for the COVID-19-like epidemics tested here (electronic supplementary material, table S2).
To examine the impacts of masking at the individual level, for each simulated population (under a specific masking setting), we compare the (re)infection rates among mask wearers with non-wearers.Across a wide range of settings tested, the model estimates substantial lower risk of reinfection for mask wearers (figure 4a).For instance, the model estimates an approximately 50% lower risk of having greater than or equal to 3 infections for mask wearers over the 5-year study period, in the main analysis assuming an R 0 of 2.6.In addition, stratifying by mask efficacy, the estimated protection against infection increases with mask efficacy, and this dose-response is evident across all age groups (figure 4b) despite different contact patterns (figure 3a).

Discussion
Despite the high efficacy of masks in reducing particle inhalation, reported effectiveness of masking has been less than ideal.This large discrepancy has puzzled researchers and prompted doubt on the validity of masking as a public health measure.Here, our model not only generates consistent results reconciling the efficacy-effectiveness discrepancy, but also helps elucidate the underlying drivers.Model results identify factors that can lower effectiveness and highlight a key 'leakage' source via household 31 31 30 30 29 28 24 Figure 4.Estimated effectiveness of masking at the individual level (main analysis assuming an R 0 of 2.6 and no masking at home).(a) Estimated per cent reduction of repeated infections over the 5-year study period (see panel subtitles), comparing mask wearers with non-wearers in the same simulated population.Numbers in the cells show the median effectiveness estimates (computed per equation (4.3)) for the corresponding masking rate (y-axis) and timing of masking (x-axis, percentages show population infection prevalence when masking is enacted).Bold fonts indicate the 2.5th% of the estimate is above 0 (i.e.significant at a = 0.05 level); note the 95% intervals here are probably narrower than measures in RCTs as the model assumes all infections are identified for the entire population.Negative values (i.e.increases; see grey cells) are possible due to model stochasticity (N.B., none are significant at a = 0.05 level).(b) Estimated per cent reduction of total infections over the 5-year study period, by age group and mask efficacy.Boxplots show estimated median (middle bar), interquartile range (box edges) and 95% intervals (whiskers).
royalsocietypublishing.org/journal/rsif J. R. Soc.Interface 21: 20230666 transmission.Further, we show the interaction of nonlinear infectious disease transmission dynamics, control measures (e.g.masking and contact reduction/social distancing), and resulting potential measurement biases.As such, the model findings here help provide a more complete picture of masking at both the individual and population levels.
For individuals, modelling results here clearly demonstrate that consistent masking can reduce the risk of initial infection and, over time, reduce the frequency of repeated infection and hence the number of reinfections.For SARS-CoV-2/COVID-19, there remains large uncertainty in the health impact of reinfections.While population-level data have suggested a likely lower severity for reinfections due to prior immune protection (e.g.data from South Africa [17]), individual (re)infection severity may vary vastly.For some individuals, there may be a substantial risk of long-lasting post-acute sequelae following SARS-CoV-2 infection (i.e.long COVID [18]) and the risk of long COVID may further increase for subsequent reinfections.In the light of the unknown long-term health impact of repeated infections, continued precautionary measures (e.g.masking) that reduce repeated infections may be advisable given these potential health impacts.
At the population level, model results here also demonstrate that masking can provide direct protection to mask wearers as well as indirect protection to non-wearers, collectively reducing epidemic intensity.In addition, as for other intervention measures, timely implementation of masking is crucial to ensure its effectiveness.However, consistent with the RCTs, model results here also indicate the overall effectiveness is probably often lower than the reported efficacy measured in laboratory settings.As demonstrated in our counterfactual modelling and testing of a wide range of masking settings, the ultimate effectiveness is limited by multiple factors ( particularly, household transmission).It is thus important to recognize this less-than-ideal effectiveness, in order to more effectively devise public health measures (e.g.planning for additional preventive measures).
To focus on the general impact of masking and influencing factors, here we sidestepped several processes shaping SARS-CoV-2/ COVID-19 dynamics; these include the emergence and circulation of different variants throughout the pandemic, differential protection of vaccines against different variants, and potential changes in disease severity for subsequent repeated infections.In particular, for more infectious variants (higher R 0 ), it may be more challenging to reduce transmission via masking, as preliminarily shown in our sensitivity analyses (electronic supplementary material, figures S1 and S3).However, the model developed here provides a general framework to study these questions and beyond.
Despite these limitations, the quantitative modelling developed here helps tease apart the various factors affecting masking effectiveness in the real world.Overall, our results suggest it is important for individuals to use masks during an epidemic, and for policymakers to recognize the non-ideal effectiveness of masking when devising public health measures.

Data sources and processing
Here, we used contact survey data from the POLYMOD study [15,19] to parametrize individual contact patterns (figure 2a), and mobility data (figure 1c) and vaccination data (electronic supplementary material, figure S6) from New York City (NYC) to model social distancing and vaccination during the pandemic.The POLYMOD study recorded contact data for 7290 participants (97 904 contacts) in eight European countries (e.g. for each contact, household versus non-household, physical versus non-physical and duration) [15,19].To compute the individual contact duration, we pooled data of all participants, weighted non-physical contact (multiplicatively by a factor 0.2 in the main analysis and 0.1, 0.5 or 1 in the sensitivity analysis) to account for likely lower per-contact infection risk than physical contact.We then summed the weighted contact durations by contact type (household versus non-household) for each individual.Note that an individual's total contact duration per day could exceed one day because simultaneous contacts are counted separately and additively combined.Given reported differences by age group, we further stratified by age group and fit the age-grouped, type-specific contact durations to a probability distribution to formulate the agent-based model.Briefly, for each age group, we fitted the data to six distributions (i.e.uniform, exponential, normal, lognormal, Weibull or gamma) and selected the best-fit distribution based on the Bayesian information criterion [20].The specific distributions used are shown in electronic supplementary material, table S3.
During the COVID-19 pandemic, many places implemented social distancing policies to reduce population contact rates.To represent these dynamics, as in our previous work [16,21,22], we used mobility data derived from Google Community Mobility Reports [23] to model contact reductions.Here, we used data for NYC during March 2020-June 2021 to formulate the model.From July 2021 onward, we assumed there was no social distancing and increased the mobility linearly to the pre-COVID level over a three-month period and remained at that level for the remainder of the study period (figure 1c).
In addition, vaccines became available in late December 2020 and mass vaccination was rolled out subsequently.To represent the impact of vaccinations, we used data from NYC (four doses, 14 December 2020-9 March 2023, at the time of this study; electronic supplementary material, figure S6) [24] to formulate the model.The first and second doses were administered three to four weeks apart as the primary series and the third and fourth doses were administered more than 6-12 months later.As there is probably little overlap of immune protection from the preceding vaccination series/dose, we additively combined vaccinations of the primary series, third and fourth doses to formulate the daily vaccination rate in the model.Assuming similar vaccination uptakes in the future, we used the mean for the same day of the year over 1 March 2022-9 March 2023 for those under 25 (roughly a 1-year period) and 1 September 2021-9 March 2023 for those aged 25 years and over (longer period given the earlier vaccination rollout for older ages) to formulate vaccination rates for days without data.For simplicity, we assumed a 90% vaccination effectiveness for all doses.

The agent-based model
The agent-based model simulates individual infection and disease history, depending on age, contact pattern, individual preventive measures (e.g.masking behaviour), prior infection and immune protection, vaccinations and immune protection, and immunity waning.In addition, at the population level, the model incorporates infection seasonality, mobility data to represent social distancing and reductions in contact during the pandemic, and vaccination data to represent mass vaccination during the study period.We describe each component below.

Figure 1 .
Figure 1.Example simulated epidemic infection rates.For these simulations, the model assumes an R 0 of 2.6 and no mask wearing at home.Coloured lines (median) and shaded areas (interquartile ranges) show daily infection rates for simulated populations with different masking rates (see legend) and timing of masking: (a) mask wearers do so outside the home at all times or (b) only when infection prevalence is greater than 1%.For context, (c) shows changes in daily population mobility (representing reductions in contact rate due to social distancing rules) and seasonal infection risk (representing effects due to environmental conditions), used in the simulations.Mobility data for New York City during March 2020-June 2021 (shown here) were used directly as model inputs.From July 2021 onward, we assumed there was no social distancing and increased the mobility linearly to the pre-COVID level over a three-month period and then remained at that level (i.e.set to 1).Trends in seasonal infection risk were represented using a sinusoidal function (see details in Methods).
cent reduction in total infections: tallies from March 2020 through this season 12 in individuals ever-infected: tallies from March 2020 through this season 13 cent reduction in total infections: tallies for this season

Figure 2 .
Figure 2.Estimated effectiveness of masking in reducing infection rates at the population level (main analysis assuming an R 0 of 2.6 and no masking at home).Heatmaps show (a) per cent reduction of cumulative total infections from March 2020 through different seasons (see panel subtitles, same below), (b) per cent reduction of the cumulative fraction of population ever-infected, i.e. regardless of number of (re)infections, and (c) per cent reduction of total infections tallied for each season.Here, a season runs 1 year from 1 July to 30 June of the following year.Numbers in the cells show median effectiveness estimates (computed per equation (4.2)) for the corresponding masking rate ( y-axis) and timing of masking (x-axis, percentages show population infection prevalence when masking is enacted).Bold fonts indicate the 2.5th % of the estimate is above 0 (i.e.significant at a = 0.05 level); note the 95% intervals here are probably narrower than measures in RCTs as the model assumes all infections are identified for the entire population.Negative values (i.e.increases; see grey cells) are possible due to model stochasticity (N.B., none are significant at a = 0.05 level).The black lines delineate cells with a median effectiveness estimate less than 25%.
distributions: household (HH) versus non-household (nonHH) contacts (b) assume people use masks at home: per cent reduction in total infections, tallies from March 2020 through this season age

Figure 3 .
Figure 3.The impact of household transmission.For these simulations, the model assumes an R 0 of 2.6 (i.e. the same as in the main analysis).(a) Example distributions of contact duration, formulated based on the POLYMOD survey (see details in Methods).(b) Counterfactual estimates of masking effectiveness in reducing total number of infections, assuming people universally wear masks at home.Numbers in the cells show the median effectiveness estimates (computed per equation (4.2)) for the corresponding masking rate ( y-axis) and timing of masking (x-axis, percentages show population infection prevalence when masking is enacted).Bold fonts indicate the 2.5th% of the estimate is above 0 (i.e.significant at a = 0.05 level); note the 95% intervals here are probably narrower than measures in RCTs as the model assumes all infections are identified for the entire population.Negative values (i.e.increases; see grey cells) are possible due to model stochasticity (N.B., none are significant at a = 0.05 level).The black lines delineate cells with a median effectiveness estimate less than 25%.
all age groups and all mask efficacy groups combined, per cent reduction in per cent reduction in total infections: by age group, by mask efficacy (example: 50% people use mask when prevalence is > 0.1%)